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(3JT) Abstract 

One-dimensional gapped spin chains with SU(N) symmetry are known to possess N different 
topological phases. In this paper, we introduce a non-local string order parameter which charac- 
terizes each of these N phases unambiguously. Numerics confirm that our order parameter may be 
used to extract a quantized topological invariant given a non-degenerate gapped ground state wave 
function. Discontinuous jumps in the discrete topological order that arise when varying physical 

I couplings in the Hamiltonian may be used to detect quantum phase transitions. 
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1 Introduction 



For a long time, the classification of distinct phases of matter was synonymous with Landau's theory 
of symmetry breaking. The latter applies to systems where the symmetry of the Hamiltonian is spon- 
taneously broken in the ground state, such as ferromagnets at sufficiently low temperature. Different 
phases can be distinguished in terms of order parameters which correspond to expectation values of 
local observables, e.g. the magnetization. Another hallmark of the theory is the existence of massless 
Goldstone modes if continuous symmetries are broken. 

Landau's paradigm was challenged with the advent of gapped physical systems in which distinct 
phases exist even though the ground state (or the ground states) preserves the same symmetries as 
the Hamiltonian. Such systems are characterized by topological order, a term that was originally 
coined for fractional quantum Hall systems, see p]. Since topological orders are usually related to 
discrete invariants, they enjoy protection against continuous deformations of the system. It might 
happen that protection against deformations is not solely due to topology but that it rather appears 
in conjunction with a symmetry which has to be preserved. In this case one speaks about symmetry 
protected topological orders. 

The Haldane phase of SU(2) invariant antiferromagnetic spin chains based on S = 1 is one of 
the first and best understood examples of a non-trivial symmetry protected topological phase. For 
the interactions it is custom to choose a specific deformation of the Heisenberg Hamiltonian as a 
representative. The resulting model is commonly referred to as the AKLT chain [2, 3j. While the exact 
ground state for the Heisenberg Hamiltonian is not known, the AKLT chain provides a convenient 
laboratory which allows to establish several important properties of the Haldane phase with full 
mathematical rigor. In particular, the ground states of the AKLT Hamiltonian are known explicitly, 
both for periodic and for open boundary conditions. Moreover, it could be proven that the chain has 
a mass gap and that ground state correlation functions of local observables decay exponentially. 

The evidence for the identification of the Haldane phase as a non-trivial topological phase of matter 
can be summarized as follows. First of all, open boundary conditions imply the existence of massless 
edge modes. The system thus exhibits a bulk-boundary correspondence which is widely regarded as 
a typical characteristic of non-trivial topological phases. It was later understood that the topological 
nature is due to symmetry fractionalization which allows the edge modes to carry a discrete ^-valued 
topological quantum number [H [5] . Secondly, all these features can also be observed in a characteristic 
entanglement spectrum [6] which provides a virtual realization of edges even in the presence of periodic 
boundary conditions. Finally, and most importantly for our present paper, there exists a non-local 
string order parameter [7], sensitive to a kind of diluted anti- ferromagnetic order, which allows to 
distinguish the topologically trivial from the topologically non-trivial phase. 

Various extensions of the AKLT setup to higher rank groups and supersymmetric systems have 
been considered, see e.g. [3ll8ll9l [T0l lllj. Other generalizations include g-deformations of the symmetry 
group which can be used to describe anisotropic spin chains [12\ 1X3] . In all these examples the matrix 
product (or valence bond) state formalism plays a crucial role [TU [151 US QZ]- The latter has also 
proven extremely useful in connection with the classification of symmetry protected topological phases 
in general one-dimensional spin systems [H 118] . Indeed, by now it is well known that topological 
phases can be distinguished based on the properties of (virtual) boundary modes that arise when 
the system is considered with open boundary conditions or when parts of the system are traced out. 
Matrix product states are relevant in this context since their boundary and entanglement properties 
are almost trivial to access. Also, there is a natural way to associate a so-called parent Hamiltonian 
to each matrix product state which, in turn, is realized as the ground state of the former. 

The classification results just mentioned yield the number of potential topological phases and 
an explicit way of constructing a representative Hamiltonian for each of them. However, given an 
arbitrary physical system, i.e. a Hilbert space, a Hamiltonian and a symmetry, no universal recipe 
how to recover its topological class is known at present. Since all topological properties are encoded 
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in the ground state wave functions, this is first of all due to the lack of knowledge of the latter. But 
even if the ground states are known exactly or approximately through a numerical calculation, the 
definition of a quantity which can be calculated efficiently and which can discriminate between all 
different topological phases is still an open problem. The degeneracy of massless edge modes might 
serve as a first indication but it still leaves ambiguities |19|. I20] I2T] . Even access to the full entanglement 
spectrum (including the energy and all additional quantum numbers) might not be sufficient as long 
as the contributions from the two edges cannot clearly be separated from each other. For this reason, 
the most promising route to a full characterization of topological phases seems to be the definition of 
suitable non-local order parameters. Important progress in this direction has recently been achieved 
in [221 ISO] (see also [23] ). It should be noted, however, that both papers focused on very specific, 
mostly abelian, discrete symmetries. An application of similar ideas to continuous symmetries and 
generic non-abelian groups remains to be worked out. 

In the present paper we will follow an alternative route and use it for the characterization of anti- 
ferromagnetic spin chains with SU(N) symmetry^ As has been shown in |21j . there are ./V distinct 
topological phases which can be realized in such chains. These N phases correspond to the N different 
ways, the center Ztv of the group SU(N) can be realized on possible boundary spins. Just as in the 
SU(2) AKLT chain before, the situation can be understood as a fractionalization of the physical 
symmetry PSU(N) = SU(N)/7*n in a setup with open boundary conditions. Our main result is an 
explicit expression for a string order parameter for SU(N) spin chains which can easily be evaluated 



once the ground state is known, see Eq. (3.1). In contrast to earlier approaches it is essential that 
our string order parameter is a matrix valued quantity. Instead of extracting the information about 
the topological phase from the absolute value of the matrix entries we will rather infer it from relative 
complex phases between off-diagonal matrix elements. It will be proven that the order parameter 
defined in this way is quantized and that it is sensitive to the representation class of boundary spins 
with respect to the action of Z^r. The string order parameter thus allows to extract a discrete 
topological invariant which permits to discriminate all N distinct phases of PSU(N) spin chains. It 
is important to note that the topological invariant will only change when the system undergoes a 
discontinuity. For this reason it may be used as a good (numerical) measure for the identification of 
quantum phase transitions. 

In order to check the validity and applicability of our analytical results we study the phase tran- 
sition between two topologically non-trivial phases of an SU(3) spin chain. Each of the two phases 
exhibits a subtle breaking of inversion symmetry through the spontaneous occurrence of boundary 
modes. For this reason the Hamiltonian cannot be written as a polynomial in the invariant scalar 
product Si ■ S% but rather requires the use of higher order Casimir operators. To our knowledge this 
is the first time that such Casimir operators are employed systematically in the formulation of spin 
chains. We then continue with a numerical investigation of mass gap and topological order using 
DMRG. The quantization of the topological order and its discontinuity at the phase transition, see 
Figure [6| provide a clear confirmation of our analytical predictions. 

Even though spin chains based on higher rank groups like SU (N) are unlikely to be found in real 
materials, there is a chance that the corresponding Hamiltonians can be engineered artificially using 
ultracold atoms in optical lattices |24| 125]. Also, special points in the moduli space of spin chains and 
spin ladders might exhibit an enhanced symmetry. This for instance happens for SU(2) spin chains 
which are known to possess an SU(3) symmetric point for a certain choice of the couplings [26]. It 
should be noted that string order parameters have also been suggested for other systems, e.g. ID 
Haldane Bose insulators [2Z] . Since the latter has been observed in experimental measurements |28| it 
seems natural that a similar experimental verification should be possible for SU(N) spin chains and 
the string order parameter obtained from Eq. (|3.1l). 



The paper is organized as follows. In Section [2] we start with a concise definition of the physical 

lr To be precise, in most cases of physical interest the actual symmetry is PSU(N) since we will be dealing with a 
special subset of representations of SU(N), see our discussion below. 
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setup under consideration and we introduce a few of the concepts that turned out to be useful in the 
classification of topological phases: Matrix product states and projective representations. Afterwards 
we provide a thorough discussion of the representation theory of su(N) and review the origin of the N 
distinct phases of SU(N) spin chains. Section [3] contains the main result of our paper. We introduce a 
string order parameter and evaluate it in the thermodynamic limit. In a series of arguments we show 
that the string order parameter includes discrete topological information and we identify the latter 
with the parameter specifying the topological phase of the spin chain. Finally, Section [4] is devoted 
to the numerical study of a family of SU{3) symmetric spin chains which interpolates between two 
topologically non-trivial phases. The toy model provides a clear confirmation of our analytical results. 
Some more technical parts of the proofs and a brief introduction into Casimir operators of su(3) have 
been moved to the Appendices. 

2 Preliminaries 

In this Section we introduce the notation and the structures that are used in the main part of our 
text. We start with a description of the physical setup and a brief outline of the matrix product state 
formalism. The latter is used to motivate the existence of N different phases of SU(N) spin chains. 
We then review some essential aspects concerning the representation theory of the Lie algebra su(N). 

2.1 Physical setup 

Throughout this paper we are considering spin chains which are characterized by the following data. 
The spins reside at sites fcona circular chain with periodic boundary conditions, the index running 
over the set k = 1, . . . , L. It will be assumed that the length of the chain is large but finite. The spins 
are described by operators Sk which take values in the Lie algebra su(N) and which act on on-site 
Hilbert spaces H^- The total Hilbert space % = ®A, ^fc * s ^ ne product of all on-site Hilbert spaces. 
For simplicity we will assume that all Hilbert spaces %k are irreducible representations of su(N) since 
otherwise the system would admit a more natural interpretation as a spin ladder instead of a spin 
chain. Finally, the dynamics of the system is described by a local Hamiltonian H which commutes 
with the total spin S = ^fc^fc- ^ can thus be written in terms of Casimir operators of su(N). 
The simplest Hamiltonians can be expressed as a function of Sk ■ Si (corresponding to the quadratic 
Casimir) where the dot denotes an su(N) invariant scalar product. More complicated Hamiltonians, 
e.g. involving many-body interactions or breaking the permutation symmetry between the two sites, 
can be defined using higher order Casimir operators. An example of this type will be discussed in 
Section |421 

Actually, the precise form of the Hamiltonian is not particularly important for the purpose of 
this paper since we will almost exclusively be concerned with properties of states. To be precise, 
our attention rests on the ground state \<p) of the system which will always be assumed to be a non- 
degenerate finitely correlated state [HI [T51 US] (non-degenerate at least in a system with periodic 
boundary conditions). Moreover, there should exist a gap to the first excited state, thus implying 
exponential decay of local correlation functions. Both properties, the uniqueness and the gap, should 
persist in the thermodynamic limit. 

The simplest way to realize an anti-ferromagnetic spin chain is as follows. The on-site Hilbert spaces 
are alternating between a space V and its dual V*. The total Hilbert space is given 
and the spin dynamics is described by the translation invariant Heisenberg Hamiltonian 



with nearest neighbor interactions. The coupling constant J is assumed to be positive, thereby favoring 
anti-parallel spin alignment. For the symmetry group SU(2) and V being the 5 = 1/2 representation, 
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the Hamiltonian (2.1) arises naturally from the electronic Hubbard model at half filling. However, 
with regard to the study of topological phases, the Heisenberg model is not ideal in many respects. 
First of all, apart from the overall normalization there are no free parameters in the Hamiltonian so it 
can only serve as a representative of one physical phase. Also, besides the fact that the ground state 
is not known exactly, the absence or presence of a gap has not been fully established. The absence 
of a gap is known for certain representations V [29]. For other representations, the existence of a gap 
can be proven in the limit of "large spin" using a mapping to a u-model with a topological 0-term 
[501 Ell E2] • More recently, the question of the Haldane gap has been revisited in [HI [Ml EI] • 

In order to realize different topological phases while retaining full analytic control over the ground 
state of the system, it is useful to consider modifications of the Heisenberg Hamiltonian which are 
obtained by generalizing the AKLT construction [21 [3] . These Hamiltonians arise as "parent Hamilto- 
nians" of specific matrix product states (MPS) [21 [151 EH E] ■ Since all our considerations take place 
on the level of ground states we will refrain from giving detailed expressions for the Hamiltonians. 
The only exception is a specific family of Hamiltonians with su(3) symmetry which will be the subject 
of Section [4] and which interpolates between two different parent Hamiltonians. It will be used to 
abandon the idealized environment of MPS parent Hamiltonians and to provide a numerical check of 
our ideas in a more realistic scenario. 

2.2 Matrix product states and topological phases 

Let the vectors \ik) denote an orthonormal basis of the on-site Hilbert spaces %k- Using an iterated 
Schmidt decomposition, any state \4>) of a periodic spin chain of length L can be written as [T7] 



ti(A^ h ■■■A^ iL )\i 1 ---i L ) , (2.2) 



ll,...,t£, 



with a certain set of matrices carrying three different indices, one physical and two auxiliary ones. 
Such a state is known as a matrix product state. To be precise, one has to distinguish different types 
of MPS depending on the behavior of the system in the thermodynamic limit L — > oo. If one wishes 
to describe the ground state of a critical system, the size of the matrices A™ will grow beyond any 
limit. In our current paper we are only interested in gapped systems and hence we will assume that 
the dimension of the matrices A^ (and their nature) stabilizes for sufficiently large values of L. The 
resulting infinite volume states are known as finitely correlated states |14l IT5l H~6] . Even though we are 
eventually interested in the thermodynamic limit, an accurate description of the physics of the system 
can be obtained by working with finite but large L for this class of states. In the presence of a finite 
gap, there are exponential corrections to expectation values which quickly die away if L is sufficiently 
large. 



The structure (2.2) arises naturally if one associates two auxiliary spaces H(k,L) an d T~L(k,R) to 
each physical site k such that %*, k R ^ = H(k+i.L)- This guarantees the existence of a maximally 
entangled state \Ik) = \q){q\ £ H(k,R) ®^-(k+i,L) where \q) refers to an orthonormal basis of 7-L(k,R)- 
The matrices A^ can be regarded as linear maps from T-L(k,L) ® T~^(k,R) to Hk- The state \<j>) is the 
image of the tensor product |J) = 1 1\) <g) • • • (g) of completely entangled pairs under the map 

21 = A^ ® • • • <g> A^K The application of the map 21 to the product of completely entangled pairs \I) 
effectively converts the tensor product into a matrix product. 

In the spin chains we are interested in, the physical Hilbert spaces Hk carry a unitary representation 
of SU(N). Moreover, the ground state \4>) should be invariant under the action of SU (^V")rj These two 
properties imply the existence of additional structures which are realized on the data of an MPS. Let 
R}- k \g) denote the representation of SU(N) on the space Hk- According to |34| . this on-site symmetry 



2 Note that invariance is implied automatically if the ground state is unique. 
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lifts to the auxiliary level as 



RW(g)-AW = D^(g)A^D^ k+1 \g)' 1 , (2.3) 

thereby promoting rl(k,L) an d H(k,R) to representations of ^[/(iV)^] In other words, the homomor- 
phisms should be equivariant projections from H(k,L) ®^H{k,R) to T~Lk, i.e. they should commute 
with the action of SU(N). 



In fact, a careful inspection of relation (2.3) shows that the physical Hilbert space Hk and the 
associated auxiliary spaces 'H(k,L/R) enter the discussion on a different footing. To understand this 
statement, let us for a moment assume that = D^ k+1 ~^ and that the auxiliary spaces form an 



irreducible representation of SU(N). In view of Schur's Lemma, the right hand side of Eq. (2.3) - 
and hence also the left hand side - is invariant in this case if g is chosen to be in the center Z^v of 
the symmetry group SU(N). In other words, descents to a linear representation of the quotient 
group PSU(N) = SU{N)/'L^ while no such requirement exists for the matrix D^. The latter only 
needs to implement a projective representation of PSU(N), 

D( 9l )D(g 2 ) = oj( gi ,g 2 )D( 9l g 2 ) with 9l ,g 2 ePSU(N) and u( 9l , g 2 ) 6 17(1) , (2.4) 

i.e. a representation up to phase factors. It is known that the projective representations of PSU(N) 
fall into N different classes when considered modulo obvious equivalences (see e.g. |21|). 

Analogous considerations apply if the assumption = fails. By choosing suitable repre- 

sentations of SU(N) on the auxiliary spaces, one can realize any symmetry group SU(N)/F on the 
physical Hilbert spaces Hk, where T C Zjv is an arbitrary subgroup of the center of SU(N). It can 
be shown that the group SU(N)/T has |T| distinct classes of projective representations [2T] . 

In a series of papers [U0H8], the projective class of the representation of the physical symmetry 
on the auxiliary spaces has been identified as a topological invariant of ID gapped spin chains. In other 
words, the projective class arising in the MPS representation of the respective ground states remains 
invariant upon deformation of the Hamiltonian. For the symmetry group PSU(N) the previous 
argument predicts exactly N distinct topological phases. For a general treatise on ID spin systems 
with continuous on-site symmetries we refer the interested reader to [21]. 

The different topological phases of a spin chain with a given symmetry can all be realized explicitly 



by defining suitable parent Hamiltonians. More precisely, for each MPS \4>) of the form (2.2) there 
exists a local Hamiltonian with the following two properties [17J: The state \4>) is the unique ground 
state of the Hamiltonian and there exists a gap. When considered with open boundary conditions, 
this construction will lead to gapless edge modes which transform according to the projective repre- 
sentations Bl = H(i,L) an d Br = %(i,,R)- Even though the energy of boundary states will receive 
corrections and the degeneracy with the ground state might get lost upon deformation of the Hamilto- 
nian, they will remain stable until the mass gap closes in the bulk. Intuitively, the correlation length 
will diverge at the phase transition, thus allowing the two boundaries of the spin chain to interact 
with each other and to disappear. 



2.3 The Lie algebra su(N) and its representations 

For a more detailed discussion of SU (N) spin chains and a concise formulation of our result we need 
to review the representation theory of the Lie algebra su(N) (see e.g. [351 E2])- The latter is the Lie 
algebra q of traceless N x N matrices and it is generated (as a vector space) by the matrices E ah 
with a / b and by H a = E aa - E a+1 ' a+1 . Here E ab denotes the elementary matrix (E ab ) cd = 5 ac 5 bd 
with a single non-zero entry in row a and column b. The diagonal matrices H a generate the Cartan 
subalgebra f) of su(N). The other generators E ab are called positive or negative roots, depending on 

3 Here we implicitly assumed that H(k,R) an d H(k+i,L) are dual not only on the level of vector spaces but also on the 
level of representations. This automatically implies the SU(N) invariance of the completely entangled pairs |7fe). 
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whether a < b or a > b. As a consequence, the Lie algebra su(N) admits a triangular decomposition 
= 0+ © f) © 5- i n t° positive roots g+, negative roots jj_ and the Cartan subalgebra f). As a Lie 
algebra, su(n) is generated by the positive and negative simple roots E ab with \a — b\ = 1. 

All finite dimensional representations V of su(N) are so-called weight representations in which all 
generators H a are represented by diagonal matrices py(H a ). By abuse of notation we will simply omit 
to write the map py in case it is clear that we are acting on a representation. If p G f)* one says that 
a vector v £ V has weight p provided that 

H a v = p{H a )v =: p a v . (2.5) 

The different eigenvalues p a can be assembled into a tuple p = (p\, . . . , pn-i) of Dynkin labels and 
should be regarded as physical charges characterizing the state v. A convenient basis for the space 
f}* is given by the fundamental weights u a which are dual to the Cartan generators in the sense that 
uj a (H b ) = 8\. A weight can thus also be written as p = ^ a /^ a w a . Any representation space V can be 
split into distinct eigenspaces with regard to the action of the generators H a . This leads to the weight 
space decomposition 

V = V„ . (2.6) 

In a finite dimensional representation all weights p are necessarily integral, i.e. p a E Z. The set of all 
weights forms the weight lattice P which is an abelian group under addition. 

Let us now turn our attention to finite dimensional irreducible representations. As is well known, 
the latter are labeled by weights A whose Dynkin labels A a are all non-negative integers. Such weights 
are called dominant. The set of dominant weights, denoted by P + , defines the fundamental Weyl 
chamber of the weight lattice P. Within an irreducible representation A, the different weights are 
all related by the application of roots a. The latter should be thought of as the charges of the root 
generators E ab (for a ^ b) with respect to the Cartan generators H a . Phrased differently, for each 
weight p in the representation A one has A — p E Q where Q is the root lattice which is generated by 
the (finite) set of roots a. 

A distinguished role is played by the adjoint representation in which the Lie algebra is represented 
on itself (regarded as a vector space) by means of the adjoint map X \— > &dx = [X, ■]. The non-zero 
weights of the adjoint representation are precisely the roots a. The — 1 simple roots have weights 
a a which are just the rows of the su(N) Cartan matrix A a j, = 25^ — <5i a _ 6 n . For our purposes it will 
be important that there exists a unique weight p = \ Yl a >o a = (1, • - - , 1) , the so-called Weyl vector, 
which has a scalar product (p,a a ) = 1 with each of the simple roots a a . The dual generator H p G f) 
is characterized by the property 

a a {HP) = 1 (2.7) 

for all a = 1, . . . , N — 1. This generator will play an important role in the definition of the string order 
parameter in Section [3} With the previous choice of simple roots one can find the following explicit 
expression for the diagonal entries of the matrix H p = diag(-fff, . . . , iT^ r _ 1 ) G f), 

N + 1 

H' = -£--a . (2.8) 
Indeed, one can easily check that this defines the unique traceless diagonal matrix with Ha—P[P +1 = 1, 



as required by Eq. (2.7). 

The final ingredient that will be needed below is the Weyl group of su(N). The Weyl group can be 
regarded as the symmetry of the root system. It consists of rotations and reflections which leave the 
set of roots invariant and is thus a subgroup of the orthogonal group in — 1 dimensions. For SU(N) 
the Weyl group is isomorphic to the symmetric group Sjy- Under the action of the Weyl group, the 
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weight lattice P may be decomposed into orbits. In our considerations below it will be crucial that 
each of these orbits has at least one element in the fundamental Weyl chamber of dominant weights 
P + . In other words, for each weight fj, G P one can find an element S G W such that fjf = S/i is in 
the fundamental Weyl chamber, i.e. // G P + . Note that the element S need not be unique. 

In order to derive the explicit action of the Weyl group on a weight it is convenient to switch to an 
alternative set of elements €4 (i = 1, . . . , N) which span the dual fj* of the Cartan algebra of su(N). 
Given any H G f) they are defined by ti(H) = Ha. Since su(N) matrices are traceless, these vectors 
satisfy the constraint Y^i e i = which leads to a slight redundancy when weights are expressed in 
terms of the e%. However, this disadvantage is compensated by the simple transformation behavior 
under the action of the Weyl group W which, for su(N), is isomorphic to the symmetric group Sn- 
Given any permutation a G Sn and the associated Weyl group element S a G W one simply has 

Sa(ei) = e u{i) . (2.9) 

Given this formula, we can deduce the Weyl group action on any weight \x = Yli c i e i- Note that the 
labels Ci are, a priori, only defined up to a simultaneous shift. We can nevertheless arrive at a unique 
description by imposing the "gauge" Yli c % = 0> an d the latter will be assumed from now on. With 
these conventions the new labels are related to the standard Dynkin labels \i a as 



E 

0=1 



JV ~ 1 AT 

a \ — v iV — a . 



This relation can be derived using the explicit form of the roots in terms of Dynkin labels, compare 
the expression for the Cartan matrix above. 

A second reason for using the epsilon basis is that €i(H a ) and ei(H p ) can be easily calculated. The 
former evaluates to €i(H a ) = 5j )0 — <5i >a +i> while the latter is ei(H p ) = —i. Note that a shift in the 
index i translates directly to a shift in ei(H p ). More precisely, let a m G Sat be the cyclic permutation 
defined by a m {i) = i + m (modulo N). Then 

e am{i) (H p ) = ei {H p )-m + Ne{i + m-N) . (2.11) 

Here, 9 is the Heaviside step function with 6(0) = 0. 

2.4 Classes of representations 

As was discussed in detail in our previous article [21], the representations A of su(N) (and hence 
of SU(N)) fall into different classes which can be interpreted as elements of the group P/Q, the 
quotient of the weight lattice P by the root lattice Q. In terms of its Dynkin labels the class of the 
representation A = (Ai, . . . , Ajv_i) is defined by 

JV-l 

[A] = aA « mod N ■ ( 2 - 12 ) 

a=l 

When representations are specified using Young tableaux, the class of a representation can be expressed 



as the number of boxes modulo N [21]. Even though Eq. (2.12) was introduced for highest weights, 
it can be extended to any weight since the expression on the right hand side is invariant under the 
action of the root lattice Q. In other words, \p] = [A] for any weight fi in a representation with highest 
weight A. As is illustrated in Figure [TJ equation ( |2.12 ) divides the weight lattice of su(N) into N 
different sublattices. 

The congruence class [A] of a representation A determines whether the representation can be lifted 
to quotients of the group SU(N). More precisely, the value of [A] fixes the action of the center Zjv of 
SU(N) on the representation A. Elements of Z^r C SU(N) are multiples Q k I]\[ of the identity matrix 
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Figure 1: (Color online) Visualization of different congruence classes for SU(2) and SU(3) in terms 
of colors. The shaded blue boxes are possible representatives of topological classes. 



with Q = exp(^p) denoting the fundamental N th root of unity. In the representation A, this element 
is mapped to the complex phase p({l k Ijs[) = We conclude that representations A with [A] = 

are linear representations of PSU(N) = SU(N) j7L^ . Likewise we may ask whether a representation 
A lifts to any of the groups SU{N)/'L q where Z g C Z^r is a subgroup of the center. This is the case if 
and only if [A] = mod q (instead of using mod N) [2"T] , 

The connection to the classification of topological phases comes in since representations A of 
SU(N) with [A] ^ only define projective representations of PSU(N). If the physical Hilbert spaces 
T~Lk transform in a linear representation of PSU(N), the (virtual) boundary spins might still transform 



in a projective representation of PSU(N) as was discussed in Section 2.2 The division of SU(N) 



representations into N distinct classes which is described by Eq. (2.12) in this way reflects the division 
of spin chains into N distinct topological classes. 

Let us finally establish the connection to the physical spin chains which have been discussed in 
Section |2.2| Since we shall be dealing with PSU(N) spin chains in this paper, the physical Hilbert 
spaces Hk (which are described by a highest weight A) should all reside in the trivial class, i.e. 
[Hk] = [A] = 0. On the other hand, the auxiliary spaces 7~L(k,L/R) can reside in non-trivial classes as 
long as their total class sums up to zero, [H(k,R)] = ~ \H(k,L)\- Together with the condition [H(k,R)] = 
— [H(k+i,L)] which arises from the duality constraint %*^ k ™ = H(k+i,L) this means that the projective 
class of the left and right auxiliary spaces, respectively, is constant all along the chain. 



3 A string order parameter for SU(N) spin chains 

In this Section we introduce a non-local string order parameter for SU(N) spin chains which reduces 
to the diluted anti- ferromagnetic order of Rommelse and denNijs [7] for N = 2. Using transfer matrix 
methods we evaluate the string order parameter on matrix product ground states and show that it 
may be used to extract a quantized topological order parameter. The latter is capable of distinguishing 
between the N different phases of PSU(N) invariant spin chains. 



3.1 Definition and interpretation 

Let \(p) be the unique ground state of our spin system. We will assume that the system has a symmetry 
group PSU (N) and that \(f>) is in a definite topological phase described by a constant t £ Zn (regarded 
as an additive group). Following the reasoning of Section 2.2, the constant t will be identified with 
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the projective class [H(k,R)} of the right auxiliary representations arising in the matrix product state 
representation of \<fi). 

In what follows we shall prove that the ground state expectation value (af'j ) of the non-local string 
order operator 



of, = #fex P| 

k=i+l 



'i-i i -1 ,,,, 



N 



H b (for i < j) (3.1) 



contains all information required to reconstruct the value of t. It serves as a convenient tool for 
the measurement of the topological phase of the system, even in cases where the matrix product 
state representation of \<p) is not known or where the nature of the auxiliary spaces - regarded as 
a representation of SU(N) - is unclear. In the previous formula, H p refers to the Cartan operator 



associated with the Weyl vector p (see Section 2.3). For SU(2), expression (3.1) reduces to the string 
order Sf exp (in ^ S z )Sj introduced by Rommelse and denNijs [7]. 

In the following Section it will be proven that, in the limit \i — j\ — > oo, the dependence of the 
string order parameter (ofj) on a and b converges exponentially to 

T ab = lim (afh = dj n t{a - b) with n = exp ^ . (3.2) 

The prefactor Cij can be used as a first rough indication of whether the system resides in a topologically 
trivial phase or not. In a trivial phase we will always obtain Cij = while in a non-trivial phase 
the prefactor is expected to be non-zeroQ Up to this point, the discussion completely parallels the 
analysis of the conventional SU (2) string order. For SU(N), however, the most important information 
resides in the off-diagonal entries, the complex phases i}^ a ~ b \ Obviously, the constant t entering this 
expression is only defined modulo N. In fact, as we shall see below, it takes values in Zat, just as 
desired. It characterizes the projective class according to which (virtual) edge modes transform and 
it thereby determines the topological phase of the state \<f>). Whenever Cij ^ 0, the value of t can be 
extracted unambiguously by calculating (or measuring) two different matrix elements and taking their 
quotient. For instance one immediately finds T 21 /T 11 = Q. 1 . Let us emphasize that a transition from 
one topological phase to another enforces the prefactor Cij to vanish since otherwise the parameter t 
cannot change its value. 

In the way it was introduced, the constant t G Zjv determines the projective class of (virtual) edge 
modes with respect to the minimal quotient PSU(N) = SU(N)/Zn of SU(N). In a concrete physical 
realization it might happen that the actual symmetry group is not PSU(N) but rather a different 
quotient SU(N) /Z g where Z q C Zjy. In this case, the projective classes are described by Z g , not by 
Zat, and t has to be considered modulo q, see [2T] . 



The attentive reader may wonder why the expectation value (3.2) still depends on i and j even 
after taking the limit \i — j\ —> oo. The answer is simple: The result of the calculation depends on the 
representation spaces used at sites i and j and hence on how the limit is performed. The dependence 
will disappear if the system is translation invariant. 

3.2 Evaluation 



The proof of Eq. (3.2) will proceed in two steps. We first prove the factorization of the matrix 
( a ij} = {Jii)(Jj ' r) m the thermodynamic limit, up to exponentially small corrections. This step uses 
transfer matrix techniques and it is intimately related to the matrix product state structure of In 
a second step we use the Weyl symmetry of the weight lattice to reduce J? L to a simpler expression. 
The latter is further analyzed in a third step from which we conclude that (Jfj) depends on a as 
( Jf L ) oc £l at . The case J,- r can be dealt with analogously. 



4 While we still lack a mathematical proof, there is numerical evidence for our assertion. 
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Step 1: Factorization 

To prove the factorization of the matrix T ab we express the ground state \(f>) = 21 |I) in terms of the 



maximally entangled state see Section 2.2 The possibility to write \4>) in this form is a direct 

consequence of the fact that \4>) can be written as a matrix product state. In the next step we use the 
intertwining property 

H a k % = %(Hk,L + H% )R ) (3.3) 

which expresses the physical spin operator H% as a sum of spin operators H£ L and H% R on the two 
corresponding auxiliary sites. Using the singlet property of \I), 

H a k ,R\I) = -H a k+ i, L \I) , (3-4) 

one easily sees that the phase factors in the string order operator afj cancel out pairwise except for 
the two boundaries. We then immediately find 

{aij} ~ m ~ m ' (3 - 5) 

where the two operators J°" L and J b - R are defined by 

jIl = (hIl + h Ir)^~ HIr and J U = n "'' ( h Il • H U) • ( 3 - 6 ) 

Note that each of these operators acts locally on two auxiliary sites. However, neither of them can 
be lifted to an operator acting locally on physical sites, i.e. there is no way to commute them back 
through 21 without rebuilding the original non-local string. 

Now that we could eliminate the non-local string connecting the two sites i and j we can evaluate 
the string order parameter using standard transfer matrix techniques [14J. For that purpose, we write 

(<t>\%Ji,LJj,R\I) = tT{E^---E^Ef l E^---E^E^E^---E^) , (3.7) 

= tr(E^---E^) , (3.8) 



where 



(4 fcl W = T,( A[k %( A[k XM x \M (3-9) 

and E^ = E^\ A pictorial interpretation of the expectation value is provided in Figure^ The two 
traces can be evaluated by diagonalization of the transfer matrices E^ k \ considered as an operator 
mapping matrices on the auxiliary space 1~L(k,R) to matrices acting on the auxiliary space Q In 

the thermodynamic limit and with \i — j\ — > oo, the only contribution will come from the highest 
eigenvalue. All other contributions are suppressed exponentially due to our mass gap assumption. It 
can easily be seen that the identity matrices on the auxiliary spaces are left and right "eigenvectors" 
with eigenvalues dim^fc/ dim"H( fc i ) and dim%fc/dim%( fc ^), respectively. Indeed, due to Schur's 
Lemma we have 

£?/UV = C R l a p and 1^4^ = C L V . (3.10) 
Moreover it is obvious that 

dim(?4) = l a/3 4 fc /UV = C R iT H {KR) W = C R dim(H (k>R) ) , (3.11) 



3 ln other words: We read the operators E from right to left. 
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X 




Figure 2: Sketch of the one-site transfer matrix E^3 (left diagram) and of the expectation value (3.7) 
(right diagram). On the right hand side, periodic boundary conditions are assumed. 



and similarly for Cl- Since E^ is a completely positive map it is guaranteed that there is no greater 
eigenvalue ([37], Prop 3.6). On the other hand, due to our mass gap assumption this eigenvalue is 
non-degenerate, even in absolute value [T7j . 

In the limit of large separation, \i — j\ 3> 1, we can rewrite the desired expectation value in a 
factorized form as 



(1 H \Ef a \l n ) (^n u jEf b \l n ) 
W " dim% dml^ " ^)( J j,R) ■ ( 3 - 12 ) 



3 

We note that the result still depends on the representation spaces describing the start and the end 
point of the original string. 

Step 2: Employing Weyl symmetry 

In the second part of the derivation we focus on the a dependence of the expectation value J a := (Jfr) 
which we claim to be proportional to Q, at with t = [Huja]. The same reasoning can be used to derive 
that (Jjn) oc 6 ^W' £ )1 = il~ bt from which the main result, Eq. (3.2), follows. Here we used the chain 
of equalities [H(j, L )] = [H*u R ^] = ~[H(i,R)] = 

Since the operators Jf L contain Cartan elements only, their expectation value can be calculated 
most easily in an orthonormal basis \ac(3) of the auxiliary space Hul) ® ^(i,-R) wri ich respects the 
weight space decomposition. In such a basis the operator J? L is represented by a diagonal matrix 
with components J£g- In order to keep the notation simple we shall use the abbreviation a & fj, if |a) 
is contained in the weight space with weight fi (of Hul,) m this case). Moreover, we wish to recall that 
the matrices (^4^)*^ are SU(N) invariant projections from auxiliary space to physical space which 
can be represented as the matrix element (s\af3). From the definition of the expectation values J a in 



Eq. (3.12) we immediately conclude 



= as^Ek'W)!'-* = xkwX £ l<sM)|V ^ ' (3 ' 13) 

In the second equality, instead of summing directly over all basis vectors, we first sum over weight 
spaces followed by a sum over vectors spanning a certain weight space. We also used an obvious 
selection rule for the weights entering the Clebsch-Gordan coefficients (s\af3). The values J^n do not 
directly depend on a and (3, but only on the weight space they belong to. We may thus define 

J % ■= Ja/3 wi t h ot E /j,, f3 <E v . (3.14) 
Furthermore, it is convenient to introduce the abbreviation 

p(ji,u) : = Yl M"/ 3 )! 2 ( 3 - 15 ) 
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such that above expression can be written as 



At this point we split the sum into orbits with respect to the Weyl group. To be more precise, we 



simplify Eq. (3.16) by restricting the summation to those weights \x and v such that their sum is in 
the fundamental Weyl Chamber, fi + v £ P + . All the other terms are obtained using the action of 
the Weyl group. Since the weights at the boundary of P + are invariant under a subgroup of the 
Weyl group this leads to an overcounting which is compensated by dividing through the order of the 
stabilizer subgroup W^+ u C W. This procedure yields 

J ' = dnk E Wh\ E P(S »- S ^S M (3-17) 
ti+ueP+ 

dirn^ E ^ E J 5M,%) ( 3 - 18 ) 



1 E S^H • ("9) 



dim Mi |W„+„ 

In the second equation the Weyl invariance of P(fJ., v) is used: v) = P(Sfi, Su) for all S £ W. 
This is proven in Appendix [A| The third equation defines K® u . Since P{^l,u) is independent of a, 
we are left to show that if" (o) oc f2 a M = $7 ai . The identification of [u] with i follows since the label 
[v] is the same for all weights v appearing in the decomposition of the su(N) representation Jium- 
This is a direct consequence of the fact that the ground state \4>) was assumed to be in a well-defined 
topological phase. 

Step 3: Weyl group gymnastics 

At this point, all ingredients are set to show that oc Jl a M, where is defined by Eq. (3.19): 



Kl v = S(fi + u){H a )n- s ^ HP ^ . (3.20) 
sew 

Writing the weights as v = aei and \i + v = dj€j (with the "gauge fixing" q = dj = 0), 



respectively, and using the Weyl group action specified in Eq. (2.9) allows us to rewrite this expression 
in the form 

K w = E* E W-E^w^-E* E ^(-E^wcW 

i <r(i)=o 3 i <x(i)=<i+l 3 ' (3 21) 

= E^(^ a) -^ a+1) ) • 

i 

Let us now focus on the sum over the different permutations a which has been abbreviated by 
in the previous formula. Our goal is to remove the constraint involving the index a in the summation 
over the permutations and to convert it into an explicit dependence of the whole expression. This can 
be achieved by the following simple resummation which makes use of the cyclic permutations a a and 
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of Eq. (2.11) 



<r(i)=o V 3 J 

E fiA (-E c -Hi)^ ( i)(^)) 



<r(i)=iV 



E ^(-E c ^o)(^ p )+«E^- E Nc ^u)) 

'i)=N V j j j=N~a+l J 



(3.22) 

AT 



The last expression arises from the following simplifications on the third row. The first sum in the 
exponent of the last equation is independent of a. The second sum vanishes due to our choice of "gauge 



fixing" Y2j c j = 0- The third sum contains the information we are after. Formula (2.10) implies the 
relation Nci = — ^ a av a = — \y\ (modulo N) for any index i. Hence each term in this sum is equal to 
[u] = t. Moreover, there are exactly a of these terms in this sum. This gives the desired dependence 
of K® v oc J7 a M on the index a. The equation 

= n-(i-tf) „ pu^ lm 

dimW, |W„ + „| ^ ' V - 



follows immediately. We have thus confirmed that the string order parameter is given by Eq. (3.2) 
and that it is a suitable tool for measuring the topological phase of a state on a spin chain. Moreover, 
the previous equation also implies that the string order operator always has a vanishing expectation 
value as long as the edge modes transform according to a linear representation of PSU(N), i.e. when 
the system is in a topologically trivial phase with t = 0. 



3.3 Properties 



The string order parameter that is derived from Eq. (3.1) has a number of desired features that one 
expects for a quantity capable of measuring a topological property. First of all, the factorization (3.12) 
implies the invariance under arbitrary block renormalization between the end points in questions. From 
a mathematical perspective this is the analogue of invariance under continuous deformations or choice 
of metric. Even though the factorized expression resembles a local correlation function one should 
bear in mind that the invariance of the ground state under PSU (N) leads to a subtle entanglement 
which propagates from site to site and cannot be removed by block renormalization [4J. 



It should be emphasized that our string order parameter (3.1) gives a reliable answer about the 
precise type of the topological phase. In contrast, entanglement entropies and spectra only encode 
information about the number of massless edge modes but not (at least not directly) about their 
representation type (see e.g. [E]). Indeed, even when only considering irreducible representations 
of SU(N), the dimension is not sufficient to distinguish between a representation and its dual for 
instance. A systematic search for even more convincing examples already succeeds for SU(3): This 
group has four different 15-dimensional irreducible representations labeled by (2, 1) and (4, 0) as well 
as their conjugates. While (2, 1) and (4,0) belong to the class [1] 6 Z3, the representations (1,2) and 
(0,4) belong to the class [2] £ Z3. So, even when forgetting about the possibility to form direct sums 
of irreducible representations we recognize that the dimension of a representation alone might not be 
sufficient to specify the topological phase it is associated with. We wish to stress that the evaluation 
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Figure 3: A graphical representation of the ground state \4>i). The black lines represent the SU(3) 
invariant mappings from auxiliary space 3 (g) 3 to physical space 8. The gray lines represent singlet 
states on 3 <g> 3. 



of the string order parameter is very efficient and fault tolerant, at least once an MPS representation 
of the ground state has been found, using DMRG for instance. 

The formula we derived for the string order and its interpretation in a sense assumes an ideal 
measurement. The form of the outcome and the particular dependence of the complex phase factor 
on the label a rely on a very specific and fixed choice of basis for the Cartan generators. In a real 
physical measurement in a laboratory one will generally measure the expectation value for a linear 
combination of operators which slightly deviates from H a . A more detailed analysis of this effect, just 
as of finite size corrections, is beyond the scope of the present article. 

4 Numerical verification 

In this Section, it will be verified in a concrete physical setup that the string order parameter defined 
in Section [3] is capable of measuring the topological order of a spin chain. For this purpose we 
define a family of SU(3) invariant Hamiltonians which smoothly interpolates between two distinct 
topologically non-trivial phases. We determine the ground states as well as the mass gap numerically 
using DMRG and we study the behavior of the string order parameter and its associated topological 
order parameter t. The numerical results clearly confirm our theoretical predictions. The complex 
phase of the string order parameter is quantized and jumps at the phase transition. 

4.1 Setup and idea 

In what follows, we shall consider a family of PSU (3) invariant spin chains with open boundary 
conditions. The on-site Hilbert spaces are all chosen to be equal to the eight-dimensional adjoint 
representation of SU(3), which is described by the highest weight (1,1). Since [(1,1)] = 0, this is 
clearly a representation of PSU (3). Moreover, we add one dynamical boundary spin at each end. 
This should be regarded as a technical trick which allows to eliminate potential edge modes thereby 
ensuring the existence of a unique ground statej^] We start with a discussion of two particular states 
\4>i) and 1^2) and their associated parent Hamiltonians H\ and H2- For these two systems we have 
full analytical control over all relevant properties such as the energy gap and the topological phase. 
We then consider the family of Hamiltonians 



and determine the ground state and the string order parameter numerically as a function of c. 

The state \4>i) is a matrix product state defined as follows: As the left and right auxiliary spaces we 
choose the two distinct three-dimensional representations 3 and 3 of SU(3), with highest weight (0, 1) 
and (1, 0), respectively. The matrices A correspond to the SU(3) invariant projections A : 3 (g> 3 — > 8 
as described in Section [2] By construction, the state \4>\} resides in the non-trivial topological phase 

6 The same effect could be obtained by imposing periodic boundary conditions which, however, are more cumbersome 
to implement in DMRG. 



H(c) 




with 



c G [0, 1] 



(4.1) 
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Figure 4: Construction of the two-site 



t = [(1,0)] = 1. As is well known, the parent Hamiltonian for an open chain of this form will lead to 
massless boundary spins transforming in the representations 3 and 3, respectively. In order to absorb 
these we couple them to extra boundary spins which transform in the dual representations 3 and 3, 
respectively (see Figure [3]) . For a suitable choice of boundary Hamiltonian we thereby end up with a 
unique ground state. A state which resides in the topological class t = 1 necessarily breaks inversion 
symmetry since the representations 3 and 3 in the auxiliary space need to be treated on a different 
footing]^] Since the ground state is required to be non-degenerate, this actually provides an interesting 
challenge for the construction of a suitable two-site Hamiltonian as will be discussed below. 

The state is obtained from \<j>) by inversion. In particular, the left auxiliary space of each site is 
interchanged with the right auxiliary space. Moreover, the reflection also includes the extra boundary 
spins which leads to a small subtlety to be explained below. As should be clear from the exchange of 
auxiliary spaces, the new state \<f>2) resides on the non-trivial topological phase t = [(0, 1)] = 2. Of 
course we can also apply the inversion to the Hamiltonian Hi, resulting in a new Hamiltonian H% of 
which |i^2) is the unique ground state. 



4.2 A family of Hamiltonians 

We are now making the preceding construction more explicit, following the standard strategy of the 
AKLT construction [2j [3] . Our goal is to find concrete expressions for the Hamiltonians Hi and H2 as 



well as for the interpolating Hamiltonian H(c) defined in (4.1). This requires introducing the concept 
of Casimir operators (see also Appendix [B]) and the calculation of a few tensor products. It turns out 
that we can restrict our attention to Hamiltonians involving nearest neighbor interactions only. 
In the bulk of the chain, the two-site Hilbert space decomposes as follows, 

(1,1)® (1,1) = (0,0)e(l,l) s e(l,l) o ® (3,0)8(0,3)0(2,2) . (4.2) 

The subscripts in (1, l) s and (1, l) a refer to the symmetric and to the anti-symmetric part of the 
tensor product. Schur's Lemma implies that su{3) invariant Hamiltonians cannot change the type 
of representation. This leaves one parameter for each of the representations which occur with mul- 
tiplicity one but four parameters for the representation (1, 1) which appears with multiplicity two. 
The latter can be thought of as the entries of a 2 x 2 matrix which acts on the multiplicity space of 
the representation (1, 1). In total, there is thus an eight-dimensional space of two-body Hamiltonians 
which commute with the action of su(3). In what follows, we will express these explicitly in terms of 
invariant combinations of the spin operators Si and S2 on the two sites. 



7 Note that the asymmetry should already be incorporated in the bulk, not just by the addition of boundary spins 
which explicitly break inversion symmetry. Otherwise our considerations would not generalize to periodic boundary 
conditions. 
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Figure 5: Construction of the two-site Hamiltonian at the left and at the right boundary with extra 
boundary spins transforming in the fundamental and the anti-fundamental representation, respectively. 
Similar considerations apply if the two boundary spins are exchanged. 



The basic objects we have at our disposal are the expression Q12 = S\ ■ S% which is related to the 
quadratic Casimir (Si + S2) 2 as well as the cubic terms C112 = drstSlSfS^ and C122 = drstS^S^S^ 
which are defined using a symmetric invariant rank three tensor d rs t, see Appendix ([B]). In addition, we 
need to consider polynomials in these objects, potentially with exchanges in the order of the operators. 
One example for such an operator would be 

^ '.= d rs t d uvw Si S\ S± 5*2 S2 S2 • (4-3) 

A careful analysis shows that the eight-dimensional space of invariant operators acting on the tensor 
product (1, 1) ® (1, 1) is spanned by 

( 1 , Q12 , Ql 2 , Q12 , C s = di 2 + C 122 , C a = C112 - C 122 , , [C a , C&] > . (4.4) 



The action of some of these operators on the constituents of the tensor product (4.2) is summarized 
in Figure |U Note that C a is an operator which exchanges the symmetric and the anti-symmetric part 
of the tensor product. After some linear algebra, it turns out that a good choice for the interpolating 
two-site Hamiltonian entering (4.1) is given b}|^] 



^Buik(c) = 1 + ^Si ■ S 2 - ^(Si ■ S 2 ) 2 - . 5 2 ) 3 + (1 - 2c) 2 -C a - ^C® . (4.5) 

We note that the deformation parameter c only multiplies the term C a which explicitly breaks inversion 
symmetry. It is not obvious at all, but an explicit calculation shows that the Hamiltonian above reduces 
to a projector for c = and for c = 1 (see the table in Figure [4]). In both cases it projects onto the 
subspace generated by (3, 0) © (0, 3) © (2, 2) as well as two (different) one-dimensional subspaces in 
the two-dimensional multiplicity space of (l,l) s ffi(l,l) a . The latter single out a specific copy of (1,1) 
inside of (1, l) s © (1, l) a . In other words, the space of zero-energy states (for two sites) is given by 
(0, 0) and states in a complementary copy of (1, 1) within (1, l) s © (1, l) a for c = and c = 1. This is 
precisely the content of (1, 0) © (0, 1), i.e. the contribution of the four auxiliary sites with the singlet 
constraint imposed, thus showing that the Hamiltonians H(0) and H(l) are of AKLT-type. 

On the boundary of the chain, the previous argument needs to be modified, compare Figure [5] We 
focus on the state \4>i) first, corresponding to c = 1. Due to the presence of extra boundary spins, the 
two-site Hilbert space now decomposes as 

Left end: (1, 0) <g> (1, 1) = (1, 0) © (0, 2) © (2, 1) (4.6) 
Right end: (1,1) ©(0,1) = (0, 1) © (2, 0) © (1, 2) . (4.7) 



8 When making this specific choice we used input from the boundary (see below) in order to fix the projection in the 
two-dimensional space (1, 1). 
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Figure 6: (Color online) A phase transition between two topological phases. In the left panel the 
energy of the ground state (blue, squares) and first exited state (red, circles) is depicted as a function 
of the interpolating parameter c. The inset clearly shows a level crossing at c = 0.5. In the right panel 
the parameters t (blue, squares) and R (red, circles) are depicted. We used adapted increments of 
0.05, 0.005 and 0.0005 in the intervals [0,0.45], [0.455,0.495] and [0.4955,0.5045], respectively. 



In both cases there is a three-dimensional space of Hamiltonians which commute with su(3). For the 
Hamiltonian we choose the projector onto the orthogonal complement of the boundary spin represen- 
tation in either case. A brief calculation, see Figure [5j shows that the latter is given by 

9 1 - - 2 

^Boundary = „ ~ r (Si ' S2) ■ (4.8) 
o o 

Note that the Hamiltonian naively does not depend on whether we are considering the left or the 
right boundary since the eigenvalues of the quadratic Casimir cannot distinguish a representation 
from its dual. However, the spins are of course represented by different matrices on both sides. The 



total Hamiltonian H(c) of the system is just the sum of Eq. (4.5) for each neighboring sites and the 
weighted (with c and 1 — c) linear superposition of the four boundary contributions Eq. (4.8). Four 
contributions arise since we have two boundaries and two choices of boundary spin representations at 
our disposal. 

The last statement also brings us back to the subtlety mentioned earlier. The spin operators in the 
Hamiltonians Hi and H2 involve different matrices for the extra boundary spins. In the interpolating 
family one thus uses two different types of spin operators (one for 3, one for 3) on one and the same 
space. Of course this breaks SU(3) invariance of the full chain. A conceptually cleaner approach would 
be to use periodic boundary conditions in the numerical simulation to be discussed below. However, 
since the extra boundary spins only lead to a localized disturbance of the wave function and since we 
are testing a property of the bulk system, we do not expect significant modifications in the outcome. 
And indeed, our numerical results clearly confirm the existence of topological order. 

4.3 Evaluation of the topological order parameter and discussion 



For different values of c in Eq. ( |4.1[ ) , we have calculated the ground state and the first excited state 
using DMRG techniques |38|, I39| . We have considered a chain of length 12 and worked with an auxiliary 
space of dimension D = 30. Figure |6(a)| shows the energy of the ground state and the first excited 
state as a function of the interpolating parameter c. We clearly recognize a gap at c = and c = 1, 
verifying that the ground states of Hi and H2 are unique. The gap tends to zero at the value c = 1/2 
where the Hamiltonian restores inversion symmetry. The closure of the gap implies the existence of a 
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phase transition. The inset clearly suggests a level crossing and hence the realization of a first order 
phase transition. 

Besides energies, we also calculated the expectation value of the string order parameter (er?- ) 
numerically, for the specific sites i = 3 and j = 9. We compared the resulting matrix to the expression 

Ja^jj = -R Q t " ^ with u = exp(2vni/3) , (4.9) 

which is the theoretical prediction for the string order parameter in the limit of an infinite chain 



(see Eq. (3.2)). The numerical values of the parameters R and t have b een e stimated by minimizing 

and they are in perfect 



6(b) 



tr(dcr • dcrt), with da = (afj 1 ) - (cr? 6 ). The results are plotted in Figure 
agreement with the theory. The parameter t is quantized and restricted to the numbers 1 and 2, thus 
providing the desired label for the topological class of the system. Moreover, this parameter changes 
discontinuously at the value c = 1/2, precisely at the point where the phase transition occurs. Small 
deviations from quantization and the failure of finding R = at the phase transition are probably due 
to finite bond dimension and finite system size. Indeed, there are finite size corrections which have 



been neglected in the derivation of Eq. (3.2). These finite size effects become more important as the 
mass gap goes to zero and the correlation length increases. Of course, the transfer matrix method 
allows to compute the string order exactly, even for finite size of the system. However, our numerical 
analysis here should merely be regarded as a proof of principle. A more accurate treatment requires 
additional work. Despite our rather brute approximation we still clearly see the crossover from one 
topological phase to another. 

We would like to remark that the numerical evaluation could equally well be performed in a system 
with periodic boundary conditions or with an arbitrary choice of ground state (in case of degenerate 
edge modes). Since all ground states coincide in the bulk we would not expect any substantial 
differences compared to the results presented here. 



5 Conclusions 

In our paper, we have searched for a physical observable which allows to distinguish the N different 
topological phases of SU(N) spin chains. To achieve this goal we have proposed a non-local string 
order operator in equation ( |3.1| ) and we have shown that its expectation value provides an unambiguous 
measure for the topological phase the chain resides in. In essence, our string order parameter extracts 
the projective class of the representations according to which potential (virtual) massless boundary 
modes transform in. It should be emphasized that, in contrast to earlier studies, our string order 
parameter is matrix valued. All matrix entries are equal in absolute value and identical to zero in 
the topologically trivial phase. The information about the - quantized - topological phase of the 
chain is contained in the relative complex phases between different matrix entries. More precisely, 
the quotient of two suitably chosen matrix elements is completely sufficient in order to extract the 
quantized topological order parameter determining the topological phase. Our analytical results are 
supported by the numerical study of a family of SU (3) Hamiltonians which interpolates between two 
distinct non-trivial topological phases. Since the realization of these two phases enforces the breaking 
of inversion symmetry, the Hamiltonian employs a new construction scheme making explicit use of 
higher order Casimir operators. We find full agreement between our analytical predictions and the 
numerical results. Indeed, Figure [6] clearly exhibits the (approximate) quantization of the topological 
order parameter, with a jump precisely occurring at the phase transition. 

It remains to be clarified how our string order parameter relates to other recent proposals for the 
determination of the projective class of (virtual) edge modes from a given matrix product ground 
state [221 [20] . While generally applicable in theory, considerable efforts need to be made in order to 
adapt these formalisms to the concrete symmetry at hand, especially when the latter is described by a 
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non-abelian group. The order parameter proposed in our paper, on the other hand, has the advantage 
of being very explicit and of being directly computable, without the detour of determining an MPS 
representation for the ground state under consideration. 

Even though tentative results have been included here, we believe that SU(N) spin chains deserve 
further numerical study. First of all, our numerical investigation only covered a special family of 
SU(3) spin chains, leaving room for an extension to higher values of N. Moreover, we made use of a 
rather brute trick in order to ensure the existence of a unique ground state: We simply absorbed the 
boundary spins which arose as a consequence of the non-trivial topology in the bulk by adding explicit 
boundary degrees of freedom. This addition could only be achieved at the cost of breaking the SU(3) 
symmetry explicitly. Since the symmetry breaking terms are localized at the boundary we believe 
that the physics in the bulk and hence the validity of our conclusions is not affected. Nevertheless it 
appears desirable to repeat the analysis with periodic boundary conditions and with higher accuracy, 
i.e. with an increased bond dimension and larger system sizes. With periodic boundary conditions the 
existence of a unique ground state can be guaranteed while preserving the full symmetry. 

One important open problem in this context is the identification of the phase transitions when 
crossing the boundary between two distinct topological phases. At least in the SU(3) example ana- 
lyzed in Section |4.3[ there is clear evidence that the phase transition is of first order, in contrast to 
potential expectations. It remains to be verified whether this statement lifts to SU(N) and whether 
the Hamiltonian can be tuned to realize a second order phase transition. Since SU(N) spin chains 
have N distinct topological phases, we expect a complicated phase diagram with a large number of 
different phase transitions which might be implemented. It would be interesting to investigate whether 
each pair of mutually distinct phases is directly connected or whether they are only connected via a 
series of phase transitions each of which changes the Zjv topological order by one unit for instance. 

It is evident that systems which are invariant under continuous symmetries different than SU(N) 
should also admit a string order parameter similar to the one described in the current paper. Even 
though the groups of type SU(N) are the most interesting ones due to the large size of their center, 
it is known plj that two and three distinct non-trivial topological phases, respectively, also exist for 
the symmetry groups E% and Spin(2N) (the universal cover of SO(2N)). Just as for SU(N) a single 
expectation value will not be sufficient to distinguish between different types of topological order 
for such symmetries. In addition, an extension to certain classes of supersymmetric or anisotropic 
systems looks feasible. It should be noted, however, that the respective symmetries of these systems 
are described by supergroups or quantum groups and that a classification of topological phases is 
still missing in that context. Nevertheless, it seems likely that our formula (3.1) will be applicable in 
anisotropic spin chains with SU q (N) quantum group symmetry without modification. 

Let us finally address an interesting conceptual issue that arises in connection with our work. For 
the original SU(2) AKLT chain it is well known that the existence of a non-trivial Rommelse-DenNijs 
string order [7] is equivalent to the breaking of a discrete hidden symmetry Z2 x Z2 [101 EH 021 03] . 
This intimate relationship can be made manifest by means of a non-local transformation of the spin 
chain. It would be very interesting to investigate whether a similar relationship exists for general 
SU(N) spin chains and to analyze the symmetry breaking patterns of discrete groups that arise in 
this way when considering the full hierarchy of topological phases |21j . The relationship between string 
order and discrete hidden symmetries for higher rank groups was also discussed in |44j . 
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A Weyl group invariance of Clebsch-Gordan coefficients 



In this Appendix it will be shown that the expression P(p, v) defined in Eq. (3.15), is invariant under 
a Weyl transformation of the weights p and v. Note that the expression can be rewritten as a trace 
over three orthogonal projections: 



PM = Yl K a l*i>f = tr(n^n^n^) . (a.i) 



Recall that Hk C T~i(k,L) ®T~^(k,R)'- denotes the orthogonal projection onto this subspace. U^ u 
denotes the orthogonal projections on the weight space <g> V u C T~i(k,L) ® 7~L(k,R)- 

The Weyl group is not only the symmetry group of the root system, but it can also be defined as 
the quotient group of the normalizer of the maximal torus with the centralizer of the maximal torus: 
W = N(T) /Z(T). The maximal torus of SU (N) simply consists of all diagonal matrices with elements 
of U(l) on the diagonal and determinant 1. Just like elements in the Cartan subalgebra h*, elements 
of the maximal torus have a simple action on states v\ with a well defined weight A. For h = expH 
(H G f)) one simply obtains 

p{h)v x = expA(#H . (A.2) 

Let a : SU(N) D N(T) — > W — > Aut(ry), where Ty is the space of weights appearing in the 
representation V. Explicitly, a w p(h) = p(w~ 1 hw). Weyl invariance of P(p,,u) will follow from 

P(a w fi,a w v) = P(p,v) ■ (A. 3) 

The advantage of this approach is that since w G SU(N), the action of the Weyl group is trivial to 
implement on V\ ® Vi. Denote : SU(N) — > Vi for i G {1,2}. Using this action we aim to show that: 

Piiw^piiw)- 1 = U awf , . (A.4) 

Since if this holds, the Weyl transformed function P can be rewritten as 

p(o w p,a w v) = tr(n CTwAtCTmI/ n w n CT 

w fJL(7 w U ) 

= tr(pi 2 (w)n Ml/ pi 2 (t(;) _1 II W/ oi2(w)n MI/ pi2(^) _1 ) (A. 5) 

= tr(n^ / 9 12 (w)~ 1 n w pi 2 (w)n MI/ ) = t^n^n^n^) = p(jj,,v) , 

which shows that P(p, v) is Weyl invariant. In the second equality pi 2 = p\® p2- In the third equality 
we make use of the cyclic property of the trace to cancel the outer to maps piz^w) and p\2{w)~ l - 
In the fourth equality we make use of the fact that p\2{w) and 11% commute. We are left to check 
the validity of Eq. (A.4). Let G and let h be an element in the maximal torus. The chain of 
equalities 

p(h)p(w)v fl = p(w)p(w~ 1 hw)v fl 

= p(w~ 1 hw)p(w)v fl (A. 6) 

= a w p{h)p{w)v fl 



shows that p{w)v^ G V aw „. From this, Eq. (A.4) follows. 
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B Casimir operators of su(3) 



The Casimir elements of a Lie algebra are polynomials in its generators S r which are central, i.e. which 
commute with each of the generators. For su(3) there are two algebraically independent Casimir 
operators. One is the usual square of the spin vector S 2 . It is associated with a non-degenerate 



invariant form and can be expressed as S 



cS r S s where invariant symmetric rank two 



tensor. The second Casimir is a cubic invariant (S, S, S) = d rs tS r S s S t which can be constructed from 
a non- vanishing invariant symmetric rank three tensor d rs t- Up to normalization, all invariant tensors 
of su(3) are obtained by choosing suitable representations and by considering traces of the form 



t 



tr(S ai ■ ■ ■ S a ") 



(B.l) 



These tensors are not all independent. On the contrary, there exist algebraic relations between the 
tensors which may be used to reduce higher rank tensors to those of relatively low degree. 



For su(3) the most convenient way of finding explicit expressions for the tensors (B.l ) is to employ 
the fundamental representation in which the spin operators S r = A r /2 are proportional to the Gell- 
Mann matrices A r (see e.g. |5]). One then defines 



tr(A r A s 



2S r 



d 



rst 



1 



tr({A r ,A s }A* 



(B.2) 



By construction, k vs and d rst are manifestly symmetric. The matrices n rs and its inverse, K rs = 5 rs /2, 
serve as a metric which can be used to raise and lower indices, just as in special and in general 
relativity. The tensors which are used for the construction of the Casimir operators are K rs and 

r 7 „ „ „ 1UVW 

^rst — >^ru' h sv' h tw Uj 



Since Casimir operators commute with the action of su(3), they are represented as scalars on 
irreducible representations. With our normalization conventions, the eigenvalues of the quadratic and 
the cubic Casimir operator, 



Q = 4rz rs S r S s and C = 8d rs4 S r S s S' 

on the irreducible representation with highest weight A are given by 

2 



Qx = (A,A + 2p) = -(A? + \ 2 2 + AiA 2 + 3Ai +3A 2 ) and 



*(Ax-A 2 ) 



2 1 

-(Ai + A 2 ) 2 + -A1A2 + Ai + A 2 + 1 



(B.3) 

(B.4) 
(B.5) 



In contrast to Q\, the cubic Casimir C\ can distinguish between a representation A = (Ai, A2) and its 
dual A + = (A2,Ai). We also see that C\ vanishes on all representations which are self-dual. 
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